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ABSTRACT 


To investigate the northern Canary Current system (NCCS), results from four 
numerical experiments of increasing complexity are examined. Experiment 1, which uses 
seasonal wind forcing only, shows that, as expected, wind forcing is the key generative 
mechanism for the current, upwelling, meander, eddy, and filament structures. 
Experiments 2 and 3 have the addition of an irregular coastline. These two experiments 
show that capes are areas for enhanced upwelling, extensive filaments, maximum current 
velocities, and enhanced growth of cyclonic meanders and eddies. Also, an embayment 
like the Gulf of Cadiz is a primary region for anticyclonic meander and eddy 
development. A fourth experiment has the additional effects of thermohaline gradients 
and Mediterranean Outflow. This complex regime has features similar to NCCS 
observations, including the generation of Meddies. 
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I. INTRODUCTION 


The eastern boundary current (EEC) system off the coasts of Iberia and Morocco 
is on the northern fringe of the Canary Current coastal upwelling system, which starts at 
~10°N, south of Dakar (Parrish et ah, 1981; Bakun and Nelson 1991). Like other EBCs, 
the mean Canary Current (CC) is a broad (-1000 km), relatively slow (-10-30 cm/s), 
equatorward, yearlong surface flow extending to an average depth of - 500 m (Wooster et 
al, 1976). The portion of the CC which lies off the Iberian west coast is sometimes 
called the Portugal Current (e.g., Tomczak and Godfrey, 1994). 

Below -100-300 m depth near the coast, a poleward undercurrent is present (e.g., 
Meincke et al, 1975; Fiuza 1982). In winter the undercurrent shoals to the north, 
occasionally reaching the surface, to form a third flow component. This current, known 
as the Iberian Current (Haynes and Barton, 1990), is a narrow (-25-40 km), relatively 
weak (less than -20 cm/s) coastal, seasonal surface poleward current (Fiuza, 1982; Frouin 
et al, 1990; Haynes and Barton, 1990). Geographically, the Iberian Current is usually 
found north of Cabo da Roca (see Figure 1 for locations), but it has also been seen as far 
south as Cabo de Sao Vincente. 

Off Portugal the undercurrent can transport water from the Mediterranean Sea 
between -500 and 1500 m depth, and water from a subtropical origin in the upper layers 
(Fiuza, 1982). Off the southwestern coast of Portugal, the water from the Mediterranean 
Sea can be as shallow as -300 m depth (Ambar, 1982). 

Ambar and Howe (1979) and Ambar (1982) have documented the separation of 
the Mediterranean Outflow (MO) into two distinct cores at -7° W: a shallow core above 
-900 m depth and a deeper core at -1200 m depth. Both cores flow westward along the 
northern coast of the Gulf of Cadiz and turn poleward around Cabo de Sao Vincente. A 
third shallower (above -700 m depth) core of MO (distinct from the two main cores) has 
been traced from the vicinity of the Strait of Gibraltar northward to -38.5° N, off western 
Portugal, following the coast of Iberia (Ambar, 1982). 
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The northern Canary Current system (NCCS) is influenced predominantly by 
equatorward, upwelling favorable winds produced by the Azores High. As described in 
Nelson (1977), the Azores High, located in the northern part of the Atlantic Ocean, is a 
semi-permanent subtropical high pressure system similar in nature and behavior to the 
North Pacific Subtropical High. The center of the Azores High migrates meridionally 
with the seasons, reaching ~27° N, its southernmost extent, in March and ridging north to 
-33° N by August. This migration causes maximum wind stress values that vary 
temporally at given locations. The summer mean east-west pressure contrast between 
Portugal and the center of the Azores High is -8 mb; however, in the winter, the contrast 
weakens to -1 mb. This pressure difference results in considerably stronger northerly to 
northwesterly winds during the summer. In the winter the winds become weaker 
northwesterly to even slight southerly winds off Iberia. The shift of maximum wind 
stress also causes the upwelling favorable winds to shift from -27° N near the Canary 
Islands in January, to -43° N off Iberia by July (Fiuza, 1982). 

Recent observations have shown highly energetic, mesoscale features such as jet¬ 
like surface currents, meanders, eddies, and filaments superimposed on the broad, 
climatological mean flow of the CC and other EBCs. Satellite images of sea surface 
temperatures have shown several filaments extending off Portugal’s coast and northwest 
Spain (Fiuza and Sousa, 1989), as well as off Cape Ghir in northwest Africa (Van Camp 
et ah, 1991; Hagen et al., 1996). Other observations have also shown anticyclonic and 
cyclonic pairs of mesoscale eddies on the order of 100 km off the coast of Iberia (Fiuza, 
1984). These features have been observed during periods of predominantly equatorward, 
upwelling favorable winds and seem to be located near prominent coastline features such 
as capes. These observations provide evidence that wind forcing and coastline 
irregularities could be important mechanisms in the formation and sustainment of many 
of the mesoscale features found in EBC regions. 

Over the past few decades numerous wind forcing models of EBCs have been 
conducted, particularly for the California Current system, the classical EBC system. 
Early studies included steady wind stress (Pedlosky, 1974) and transient wind forcing 
(Philander and Yoon, 1982). The response of reduced gravity models to realistic coastal 
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winds was investigated by Carton (1984) and Carton and Philander (1984). A series of 
experiments were also conducted by McCreary et al. (1987) using a linear model with 
both transient and steady wind forcing in the California Current system. These previous 
models produced weak currents (5-10 cm/s) and no eddies or filaments developed. 

In the last decade, modeling studies have focused on the driving mechanisms of 
complex features, such as filaments and eddies observed in EEC regions. While Dceda et 
al. (1984a, b) and Haidvogel et al. (1991) studied baroclinic and barotropic instability, 
coastline irregularities, and bottom topography as possible mechanisms, Batteen et al. 
(1989), McCreary et al. (1991), Pares-Sierra et al. (1993), and Batteen (1997) studied 
wind forcing as a possible generative mechanism. Recently, Batteen (1997) used a 
multilevel primitive equation (PE) model of the California Current system from 35° N to 
47.5° N to study the effects of seasonal winds and coastline irregularities. Batteen and 
Vance (1998) studied the additional effects of thermohaline gradients. In a more recent 
study Batteen et al. (1998) studied the contribution of seasonal wind forcing, 
thermohaline gradients and irregular coastline geometry to the generation of eddies and 
filaments for the entire California Current system (22.5° N to 47.5° N). The results 
showed that wind forcing was the dominant process responsible for many of the observed 
features of the California Current system. 

In contrast to the California Current system, coastal modeling studies of the 
NCCS have been scarce. McClain et al. (1986) performed the first limited modeling 
study in the region, using ship-derived winds to produce a large negative wind stress curl 
off the northwest coast of Iberia. This resulted in opposing equatorward and poleward 
surface currents. Batteen et al. (1992) used an eddy resolving PE model with both 
uniform and variable wind stress, the latter computed from synoptic surface pressure 
analyses, to produce realistic currents as well as mesoscale eddies. Both studies, 
however, only covered a limited region of the NCCS (i.e., the northwest coast of Iberia) 
and had a straight coast in their model domains. 

While the predominantly wind-driven NCCS would be expected to produce 
similar features as in a classical EBC system, the NCCS has several distinguishing EBC 
features. It has the unique geographical feature of the relatively large embayment known 
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as the Gulf of Cadiz. It also has the MO, which spreads saltier waters outward from the 
Gulf of Cadiz, and is a source region for Meddies. 

The objective of this process-oriented study is to use a multilevel PE model to 
investigate the roles of classical as well as unique features in the NCCS. Toward this 
end, results from several numerical experiments of increasing complexity (see Table 1) 
are systematically examined. The results of a classical wind-driven EEC system along a 
straight coastline are first presented in Experiment 1. The additional effects of irregular 
coastline geometry, which includes the Gulf of Cadiz and coastline promontories, such as 
capes, are presented in Experiments 2 and 3. Lastly, in Experiment 4, the additional 
effects of thermohaline gradients and the MO are explored. 
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11. MODEL DESCRIPTION 


A. MODEL EQUATIONS 

This study uses the eddy-resolving, limited-area EEC model of Batteen (1997) 
with open northern, western, and southern boundaries. The model is multi-level, uses 
non-adiabatic primitive equations on a beta-plane, and has both baroclinic and barotropic 
velocity components. It is based on the hydrostatic, Boussinesq, and rigid lid 
approximations. Equations governing the model are as follows: 



d^u 

dz^ 

(1) 


d\ 

de 

(2) 

o 

II 


(3) 

1 

II 


(4) 

p = p,[\-a{T-T,) + p{S-S,)\ 


(5) 

dT . d^T 


(6) 

dS . d^S 


(7) 


In the above equations, t is time, T is temperature, S is salinity, p is density, and p is 
pressure. A right-handed Cartesian coordinate system (x,y,z) is used where x points 
toward shore, y alongshore, and z upward. The corresponding velocity components are 
(u,v,w). Table 2 provides a list of symbols found in the model equations, as well as 
values of constants used throughout the study. 

A space-staggered B-scheme (Arakawa and Lamb, 1977) is used for the horizontal 
finite differencing. Batteen and Han (1981) have shown this scheme to be appropriate 
when the grid spacing is approximately on the same order as, or less than, the Rossby 
radius of deformation. The horizontal grid spacing is 9 km in the east-west direction and 
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11 km in the north-south direction, while the internal Rossby radius of deformation is ~30 
km. The model uses ten vertical layers, with constant z-levels, at depths of 10, 30, 75, 
150, 250, 400, 600, 1226, 2283, and 3656 m. This vertical spacing scheme of the model 
concentrates more layers above the thermocline, which is in the dynamically active 
portion of the ocean. 

The model domain (Figure 1) encompasses the west coasts of Iberia and Morocco, 
from 30° N to 42.5° N (1408 km alongshore), and from 5° W to 17.5° W (1152 km cross¬ 
shore). The coastal boundaries of the model domain are closed, and have both the 
tangential and normal components of velocity set to zero. Bottom topography has been 
omitted to focus on the roles of climatological wind forcing, thermohaline gradients, and 
of the MO. The eastern boundary is modeled as a vertical wall. The constant depth used 
in the model is 4500 m. 

The northern, southern, and western borders are open boundaries, which use a 
modified version of the radiation boundary conditions of Camerlengo and O’Brien (1980). 
Biharmonic lateral heat and momentum diffusion is used in the model with the same 

choice of coefficients (i.e., 2.0 xlO^^ cm'^/s) as in Batteen (1997). Holland (1978) 
showed that highly scale-selective biharmonic diffusion acts predominantly on 
submesoscales, while Holland and Batteen (1986) found that baroclinic mesoscale 
processes can be damped by Laplacian lateral heat diffusion. As a result, the use of 
biharmonic lateral diffusion should allow mesoscale eddy generation via barotropic 
(horizontal shear) and/or baroclinic (vertical shear) instability mechanisms. Batteen 
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(1997) used weak (0.5 cm /s) vertical eddy viscosities and conductivities. Bottom stress 
is parameterized by a simplified quadratic drag law (Weatherly, 1972), as in Batteen 
(1997). 

The method of solution is straightforward with the rigid lid and flat bottom 
assumptions because the vertically integrated horizontal velocity is subsequently 
nondivergent. The vertical mean flow can be described by a streamfunction which can be 
predicted from the vorticity equation, while the vertical shear currents can -be predicted 
after the vertical mean flow is subtracted from the original equations. The other 



variables, i.e. temperature, density, salinity, vertical velocity, and pressure, can be 
explicitly obtained from the thermodynamic energy equation (6), equation of state (5), 
salinity equation (7), continuity equation (3), and hydrostatic equation (4), respectively 
(For more complete details on the method of solution, see Batteen, 1997). 

B. FORCING CONDITIONS AND EXPERIMENTAL DESIGN 

To explain the effects of seasonal wind, the model is forced from rest with 
climatological wind fields from a 2.5° by 2.5° grid of the European Centre for Medium 
Range Weather Forecasts (ECMWF) near-surface wind analyses (Trenberth et al., 1990). 
The monthly mean stresses based on twice daily wind analyses from 1980-1989 have 
been interpolated spatially to the 9 by 11 km model resolution and temporally to daily 
wind values. Sample wind fields used can be seen in Figure 2, which shows the annual 
migration of the Azores Subtropical High from the south in the winter (e.g.. Figure 2a) to 
its most northern extent in the summer (e.g.. Figure 2d). The atmospheric pressure 
pattern for November (not shown) and December (Figure 2a) depicts a low to the north 
and the Azores High to the south. This results in a wind divergence with weakly poleward 
winds north of -40° N and equatorward winds to the south. During January (not shown) 
and February (Figure 2b) the divergence in the wind field migrates poleward. By March 
(not shown) and April (Figure 2c) an equatorward component in the wind field is 
observed along the entire domain. The strongest equatorward winds are discernible from 
July (not shown) through August (Figure 2d). In September (not shown) and October 
(Figure 2e) the winds start to weaken throughout the domain, and divergence in the wind 
field is observed in the north by October (Figure 2e). By November (not shown) the wind 
divergence returns to -40° N. 

Table 1 summarizes the design of the four experiments. Note that in all the 
experiments the model is forced from rest with seasonal ECMWF winds. The model is 
then allowed to reach a quasi-equilibrium state by running each experiment for three 
years. 
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As Table 1 shows, in Experiment 1, to isolate the response to seasonal wind 
forcing from coastline variations, a straight coastline is used throughout the entire domain 
(e.g., see Figure 4a from Experiment 1). The initial mean stratification used are annual 
climatological temperature fields based on Levitus and Boyer (1994). The temperatures 
(°C) used for the ten levels from the surface to 4500 m are 17.5, 17.3, 16.5, 15.0, 13.7, 
12.5, 11.0, 8.52, 3.59, and 2.09, respectively, while a salinity constant of 34.7 is used for 
all levels. 

To examine the additional effects of irregular coastline geometry, two additional 
cases are run. In Experiment 2 an idealized “straightened” coastline is used with two 
prominent capes (i.e., Cabo de Sao Vincente and Cape Ghir), and an embayment 
representing the Gulf of Cadiz (e.g., see Figure 8 from Experiment 2). In Experiment 3, a 
coastline with extensive coastline irregularities (e.g., see Figure 10a from Experiment 3) 
is used to investigate the additional effects of more realistic capes on the formation and 
sustainment of the features found in the NCCS. 

Finally, in Experiment 4, the additional effects of the MO and of thermohaline 
gradients on the NCCS are explored. The salinity at each level for the MO is input as a 
constant value for five grid points at -36° N every time step with the values given in 
Table 3 (from Levitus et al, 1994). 

Seasonal temperature and salinity climatological conditions for the upper seven 
levels from Levitus and Boyer (1994) and from Levitus et al. (1994) are used to initialize 
the model, and, every third day, to force the model at the western boundary at 17.0° W. 
The seasonal temperature and salinity forcing conditions, which are initially assumed to 
be zonally homogeneous, are shown for one of the seven levels in Figure 3 for the 
northern (42.5° N) and southern (30° N) boundaries of the model domain. Since the 
lower three levels do not exhibit much horizontal variation, they are assumed to be 
constant for each level. The temperature values used for levels 8-10 are 8.52° C, 3.59° C, 
and 2.09° C. The salinity constant used for the lower three levels is 34.7 psu. 

While the temperatures to the south are warmer than those to the north above 
1000m depth, only the upper level temperature conditions (e.g.. Figure 3) show 
significant seasonal variability with a tempera,ture maximum in September and a 
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temperature minimum in February throughout the whole region. Below around 150m 
depth (not shown), both the seasonal temperature fluctuations and the temperature 
gradient weaken, as expected. In contrast, the salinity conditions above 1000m (e.g.. 
Figure 3), which in the upper levels show less (more) saline water to the north (south), 
have no significant seasonal cycle. 
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III. RESULTS FROM THE WIND-FORCED MODEL SIMULATIONS 


A. EXPERIMENT 1 - EFFECTS OF STRAIGHT COASTLINE GEOMETRY 

1. Spin-up Phase 

In response to the predominantly equatorward winds over the model domain, a 
coastal equatorward surface current develops. The current develops initially in the 
southern end of the domain in response to the stronger equatorward winds associated with 
the southerly position of the Azores High. By summer (e.g., Figure 4a) there is generally 
equatorward flow throughout the domain. In addition, a poleward undercurrent has also 
developed, initially, in the southern end of the domain and then extending along the entire 
coast (not shown). 

A cross-section of meridional velocity (e.g.. Figure 5) shows the typical vertical 
structure of the equatorward coastal jet and poleward undercurrent during the upwelling 
season. The coastal jet core remains within ~ 50 km of the coast and extends from ~ 200 
m depth near the coast to -300 - 600 m depth offshore. The core velocities range from - 
30 - 40 cm/s. Below the equatorward jet, a weaker undercurrent with a core velocity of 

- 5 cm/s at 600 m depth is also seen. The offshore extent of the undercurrent core is 
confined to - 50 km off the coast. Typical core velocities for the undercurrent range from 

- 3 to 10 cm/s. 

As the wind becomes equatorward over the entire domain and increases in 
intensity, increased stress is exerted on the ocean surface creating Ekman transport 
offshore. This results in upwelling of cooler water along the coast. In the spring the 
upwelling is confined to the south where the stronger winds exist, while in summer, it is 
present all along the coast (e.g.. Figure 4a). 

As the coastal jet and undercurrent become fully established (for example, see 
Figure 5, which shows the structure of the currents just prior to meander formation), the 
currents become unstable, due to strong vertical and horizontal shear instabilities between 
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the equatorward jet and poleward undercurrent, and form meanders. This first occurs in 
the equatorward end of the model domain where there has been more time for the strong 
equatorward winds to act on the ocean surface. As the meanders intensify, cold, 
upwelling filaments and eddies develop (e.g.. Figure 4b). These eddies are predominantly 
cold core and cyclonic, 0(100 km) in diameter, and extend ~ 50 - 150 km from the 
coast. Throughout the summer, which corresponds to the period of maximum 
equatorward wind stress (e.g.. Figure 2d), meanders, filaments, and eddies continue to 
develop equatorward of -37° N. 

In the fall, the currents begin to weaken with the slacking winds. The upwelling 
also ends, first in the poleward end of the domain and then progressively equatorward 
until by - day 330 all upwelling has ceased. By day 315 a surface coastal poleward 
current appears north of -38° N in response to the weak poleward winds present off the 
northwestern coast of Portugal in November and December. The current persists in the 
northern portion of the domain throughout the winter (e.g.. Figure 4c). South of -38° N, 
the equatorward surface current, embedded with cyclonic eddies, has taken the form of a 
meandering jet (e.g.. Figure 4c), which continues to meander downstream throughout the 
winter. 

During the second year, in the spring, with the return of the Azores High along the 
west coast, upwelling of cooler water along the coast redevelops (not shown). Throughout 
the upwelling season, cold, upwelling filaments develop inshore of cyclonic meanders of 
the equatorward jet (e.g., Figure 4d). Cold, cyclonic eddies also tend to form offshore of 
the filaments, while warm, anticyclonic eddies tend to develop downstream of the 
filaments in regions where the equatorward jet meanders anticyclonically (e.g.. Figure 4d 
shows an anticyclonic eddy centered at -13° W, -33.5° N). In the fall, as expected, the 
filaments and cold, upwelled water near the coast disappear. A surface poleward current 
reappears in the poleward end of the model domain, while to the south, the meandering, 
equatorward jet, embedded with eddies, persists throughout the fall and winter. 
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2. Quasi-equilibrium Phase 

By year three, the system has reached a quasi-equilibrium state. To see the 
seasonal structures of the features of the NCCS, the model output fields are time- 
averaged every month. 

In the spring (e.g.. Figure 6a), a relatively strong (~20-30 cm/s) coastal, 
equatorward current develops. In addition a relatively strong (~35 cm/s) equatorward jet 
is present -200-250 km off the coast. Both the equatorward current and jet meander 
downstream until -37° N where they join offshore of a cyclonic eddy to the south. Both 
the current and jet continue to meander downstream. Both cyclonic and anticyclonic 
eddies of 0(100 km) diameter are discernible both offshore and inshore of the 
meandering current and jet. 

In summer (e.g.. Figure 6b), upwelling of cooler water is present all along the 
coast, and reaches a maximum off Portugal in August. Evidence of cold offshore 
flowing, upwelling filaments south of -35° N is also found. There is both a relatively 
strong coastal equatorward flow present along the coast and an offshore equatorward jet. 
Below the coastal equatorward current there is a poleward undercurrent (Figure 7a) which 
is present all along the coast (Figure 7b). The undercurrent meanders within -100 km of 
the coast. Offshore of the undercurrent are both cyclonic and anticyclonic eddies of 
0(100 km) diameter. 

• By fall (not shown), upwelling has weakened along the entire coast. Winter 
results (e.g.. Figure 6c) show a meandering equatorward surface jet, embedded with 
eddies, all along the coast. 
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B. EXPERIMENT 2 - EFFECTS OF IDEALIZED IRREGULAR COASTLINE 
GEOMETRY 

In Experiment 1, a straight coastline was used. In Experiment 2, to examine the 
additional effects of irregular coastline geometry, the idealized “straightened” coastline 
shown in Figure 8 is used with two prominent capes and an embayment (representing the 
Gulf of Cadiz). All other parameters are the same as in Experiment 1. 

1. Spin-up Phase 

North of the Gulf of Cadiz (i.e., north of -37° N), the current and temperature 
patterns in Experiment 2 are expected to show many of the same basic features that exist 
in Experiment 1. South of —37° N, the change in coastline orientation from a north-south 
to an east-west orientation and the presence of the relatively large embayment 
representing the Gulf of Cadiz, are expected to play significant roles in the generation of 
meanders and eddies in the NCCS. 

As in the spin-up phase in Experiment 1, an equatorward surface current develops 
first in the south and then in the north as the equatorward wind forcing migrates north 
with the Azores High. While the upwelling in Experiment 1 develops fairly uniformly 
along the coast, first in the south and progressively to the north, the upwelling in 
Experiment 2 is relatively patchy and is most intense in the southern part of the domain 
(e.g.. Figure 8). The upwelling is generally weaker along the south coast of Portugal and 
the west coast of the Gulf of Cadiz. In addition, a poleward undercurrent (not shown) 
also develops, initially in the southern end of the domain and then extending along the 
entire coast. As in Experiment 1, the undercurrent has a typical core speed of -5 cm/s at 
-600 m depth (not shown). 

Once the coastal jet and the undercurrent become fully established, the currents 
become unstable and form meanders. The meanders and cold, cyclonic eddies in 
Experiment 2 are generally much more pronounced and more intense than in Experiment 
1 (e.g., compare Figures 4b and 8), and can be “anchored” to coastal promontories. Of 
particular interest is the presence of two cold cyclonic eddies which develop north and 




south of the Gulf of Cadiz at coastal promontories. The southern side of the cyclonic 
eddy (at -36.5° N, 9° W) on the south coast of Portugal and the northern side of the 
cyclonic eddy (at -34.5° N, 8.5° W) entrain the large anticyclonic meander in the Gulf of 
Cadiz with the result that an anticyclonic eddy develops in the Gulf of Cadiz. 

2. Quasi-equilibrium Phase 

To see the seasonal structures of the features of the NCCS, the model output fields 
are time-averaged every month. In the spring (e.g.. Figure 9a), a relatively strong coastal 
equatorward jet of -25 cm/s and a weaker offshore (-250 km) equatorward current of -10 
cm/s exist off the coast of Portugal. The coastal jet has a meandering pattern. North and 
south of the Gulf of Cadiz, it meanders cyclonically, while in the Gulf of Cadiz it 
meanders anticyclonically. 

By summer (e.g.. Figure 9b), upwelling is present all along the coast and in 
August reaches a maximum. Beneath the surface (not shown), a core of poleward flow 
with a speed of -5 cm/s at -600 m depth, is discernible. Since the coastal jet and 
undercurrent are fully established by summer, the currents become unstable, and form 
meanders. The coastal jet flows alongshore off Portugal until -39° N, where it develops a 
cyclonic meander. Off Cabo de Sao Vincente, the coastal jet veers offshore again to form 
a second cyclonic meander. The coastal jet then meanders anticyclonically into the Gulf 
of Cadiz. At -34° N it forms a third cyclonic meander. The largest meanders and eddies 
occur to the south in response to the stronger currents. As the strengthened currents flow 
offshore in the regions of the developed meanders, the current advects cold upwelled 
water offshore. 

In fall (not shown), upwelling has weakened along the entire coast. By winter 
(Figure 9c), a weak meandering equatorward flow is discernible. Embedded in this flow 
are numerous jets and relatively large cyclonic and anticyclonic eddies south of Cabo de 
Sao Vincente. 
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C. EXPERIMENT 3 - EFFECTS OF REALISTIC IRREGULAR COASTLINE 
GEOMETRY 

In Experiment 2, an idealized “straightened” coastline was used with two 
prominent capes (i.e., Cabo de Sao Vincente and Cape Ghir) and an embayment 
representing the Gulf of Cadiz. In Experiment 3, a more realistic coastline with extensive 
coastline irregularities (e.g.. Figure 10a) is used to investigate the additional effects of 
more realistic capes on the formation and sustainment of the features found in the NCCS. 
All other parameters are the same as in Experiments 1 and 2. 

1. Spin-up Phase 

As in the previous experiments, the current and temperature patterns in 
Experiment 3 show many of the same features. In response to predominantly 
equatorward winds, a coastal equatorward surface current develops in the southern end of 
the domain, and by day 60 extends along the entire coast (not shown). 

When the Azores High migrates north in the spring, the wind intensifies and 
transitions to equatorward flow over the entire domain. Upwelling first occurs in the 
south where the stronger winds exist (not shown). While the upwelling in Experiment 2 
was intensified on the west coast in the southern part of the domain (e.g.. Figure 8), 
upwelling in Experiment 3 is also enhanced near prominent capes to the north on the west 
coast (e.g., off Cabo da Roca in Figure 10a). Since these capes are areas where the 
alongshore component of the wind stress is at a local maximum, the coastal current 
velocity is also at a maximum, which leads to enhanced upwelling (Batteen, 1997). In 
addition, local enhanced upwelling can occur south of Cabo de Sao Vincente and off 
Cabo de Santa Maria (e.g.. Figure 10a), a feature not seen in Experiment 2 with the east- 
west orientation at -37° N. 

Below the equatorward surface current, a poleward undercurrent develops first at 
the equatorward end of the domain and then extending all along the coast (not shown). A 
cross-section of the opposing currents (Figure 11) shows the typical structure of the 
currents off the west coasts. The coastal jet extends from -200 m depth near the coast to 
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-600 m depth offshore, and has a peak core velocity of -35-40 cm/s (Figure 11a). Below 
the coastal jet along the west coasts of Portugal (e.g., Figure 11a), and Morocco (not 
shown), the undercurrent develops separate cores: one between -200-600 m centered at 
-400 m depth and one below -1000 m depth centered at -1200 m depth. Both cores have 
peak velocities of -5-10 cm/s and remain within -50 km of the coast. Off the west coast 
of the Gulf of Cadiz below -400 m depth, there is a single undercurrent core with a peak 
velocity of -5 cm/s (Figure lib). As the coastal Jet and the subsurface undercurrent 
become fully established, the flow becomes both barotropically and baroclinically 
unstable and leads to the formation of meanders and filaments. 

In Experiment 2 the majority of features had no geographical preference, except 
for the features that appeared tied to the major geographical features depicted in the 
idealized coastline (e.g., the cyclonic eddy off Cabo de Sao Vincente, the anticyclonic 
flow through the Gulf of Cadiz, and the large filament off Cape Beddouzza). In contrast. 
Experiment 3 frequently shows additional eddies and filaments “anchored” off less 
prominent coastal features (e.g., off Cabo da Roca, Cabo de Santa Maria, and the Cape 
Beddouzza/Cape Ghir area shown in Figure 10b) that occur along the realistic irregular 
coastline. In addition, in Experiment 3 upwelling occurs earlier (-2-3 weeks) and with 
greater intensity off Portugal, where the largest differences occur between the idealized 
and realistic coastlines (e.g., compare Figure 8 with Figure 10b). Experiment 3 also 
develops more meanders in the vicinity of the additional capes (which also form earlier 
than those in Experiment 2). These meanders intensify and develop into predominantly 
cyclonic cold core eddies, which in time coalesce with other eddies to form larger 
cyclonic eddies on the order of -100-300 km. Similar to Experiment 2, by the end of the 
upwelling season (not shown) the equatorward surface flow has taken the form of a 
meandering jet embedded with predominantly cyclonic eddies. 

2. Quasi-equilibrium Phase 

Similar to previous experiments, the model reaches quasi-equilibrium by year 
three, and shows that most features continue to be generated and maintained. The model 
results show that in the spring (not shown) there is a mean equatorward jet embedded 
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with numerous eddies. The jet flows close to the coast off Portugal, but leaves the coast 
as it impinges on Cabo de Roca. It then continues to meander equatorward. By late April 
the coastal jet appears off the coast of Morocco and upwelling resumes. 

In the summer (e.g.. Figure 12) upwelling is visible along the western and 
southern coasts of Iberia, as well as along the coast of Morocco. Cyclonic and 
anticyclonic eddies are discernible over the entire domain. As the equatorward jet flows 
off Portugal, it is deflected offshore by the large capes (e.g., Cabo da Roca and Cabo de 
Sao Vincente), and advects the cold, upwelled water with it, resulting in the formation of 
filaments. Maximum upwelling occurs by late summer. Beneath the surface (Figure 
13a), cores of poleward flow with speeds of ~5-7.5 cm/s at ~600 and 1200 m depth are 
discernible. An additional narrow (< 20 km wide) third core of poleward flow at -250 m 
depth with a speed of -5 cm/s is also seen off Portugal in the vicinity of -37.5® N (Figure 
13a). 

In the fall and winter, upwelling has weakened and a meandering equatorward 
current embedded with westward propagating, cyclonic eddies is discernible (not shown). 
In response to the weak poleward winds during November and December (Figure 2), a 
poleward flow develops north of -40° N (not shown). This flow, which has speeds of 
-20-25 cm/s and extends to -200 m depth, lasts throughout the winter (e.g.. Figure 13b). 
Beneath the surface cores of poleward flow at -600 and -1200 m depth with speeds of 
-10-17.5 cm/s depth are also discernible. An additional narrow third core of poleward 
flow at -250 m depth with a speed of -7.5 cm/s is also shown. 

D. EXPERIMENT 4 - ADDITIONAL EFFECTS OF THERMOHALINE 
FORCING AND MEDITERRANEAN OUTFLOW 

In Experiment 3, a coastline with extensive coastline irregularities was used to 
investigate the effects of realistic capes. In Experiment 4 the additional effects of 
thermohaline gradients and of the MO on the NCCS are explored. 
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1. Spin-up Phase 

Due to the combination of thermohaline gradients and wind forcing, different 
oceanic responses are expected depending on the season. As expected, when the 
temperature decreases poleward (Figure 3), resulting in an increase of density poleward, 
the pressure gradient due to the temperature (density) gradient establishes an onshore 
geostrophic inflow from the interior ocean. In the winter (e.g.. Figure 14a), in the 
poleward end of the model domain, the large pressure gradient and poleward wind stress 
result in onshore geostrophic and Ekman flow. As the flow approaches the coast, it turns 
and forms a poleward boundary current. In the equatorward end of the model domain, the 
smaller pressure gradient and equatorward wind stress result in weak onshore geostrophic 
flow, offshore Ekman flow, and a coastal equatorward surface current. 

During the peak upwelling season (July through September), the combination of a 
weakened pressure gradient and increased equatorward winds over the entire model 
domain lead to a strengthening of equatorward flow all along the coast and a weakening 
(strengthening) of onshore (offshore) flow. Upwelling can be seen along the coast, but 
the coldest water is found near coastal promontories (e.g., off Cape Ghir and Cape 
Beddouzza in Figure 14b). A typical cross-section (Figure 15a) of the currents shows that 
there is a surface coastal equatorward jet which extends to -180 m depth near the coast 
and -800 m depth offshore. The surface current is within -100 km of the coast and has 
core velocities of -20-25 cm/s. Below the equatorward current is a poleward 
undercurrent within -50 km of the coast which has a core velocity of -5 cm/s at -1200 m 
depth. 

In the fall, as expected, the surface poleward flow strengthens in the coastal 
poleward end of the model domain in response to both the strengthening of the pressure 
gradient and the return of poleward wind stress in the region. A typical cross-section of 
the currents in the poleward end of the model domain north of -41° N (e.g.. Figure 15b) 
shows that the coastal poleward surface current has a core velocity of -10 cm/s and is 
within -100 km of the coast. To the west of the current is the year-round equatorward 
current with a core velocity of -25-30 cm/s. 
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At depth, a relatively narrow (< 50 km wide) coastal equatorward undercurrent is 
discernible. The undercurrent is within ~30km of the coast, extends from -400-1000 m 
depth, and has a core velocity of -2 cm/s. 

The equatorward flow to the west of the poleward undercurrent below 
-900 m depth is the inshore branch of an anticyclonic eddy (Figure 14c), which is 
relatively salty (Figure 15c). It has a core velocity of - 3 cm/s at - 1200 m depth (not 
shown), and forms at - 41° N on the inshore side of the poleward undercurrent (Figure 
14c). Other anticyclonic eddies at - 1200 m depth are also discernible in Figure 14c near 
the coast of Portugal at -39.5° N, and at -37.5° N. The salty characteristics of these 
anticyclonic eddies are consistent with Meddies which have been found at -1200 m depth 
in these regions (e.g., Richardson and Tychensky, 1998). 

2. Quasi-equilibrium Phase 

Longer run times (-3 years) of the model simulation show that while many of the 
features simulated in the previous experiment are also seen in Experiment 4, there are 
some notable differences. A comparison of the temperature and currents for Experiments 
3 (Figures 12) and Experiment 4 (Figure 16) during the upwelling season show similar 
results, except for a tighter temperature gradient all along the coast in Experiment 4. 

In the winter, in the poleward end of the model domain north of -41° N, an 
equatorward undercurrent, a feature not present in previous experiments, is discernible 
(e.g.. Figure 17). This undercurrent is within -30 km of the coast, extends from -1000- 
1800 m depth, and has a core velocity of -3 cm/s. Consistent with previous experiments, 
a surface poleward current with a core velocity of -17.5 cm/s is also present. 

To clarify whether the tighter temperature gradient and equatorward undercurrent 
are due to the effect of thermohaline gradients, or to the MO, or to a combination of the 
two, a sensitivity experiment was run which was the same as Experiment 4, but without 
the effect of the MO. Model results (not shown) showed that both features were still 
present so that these features must be due to onshore geostrophic flow resulting from the 
alongshore pressure gradient. 
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Of particular note is the subsurface spreading of MO, which at -36° N, 6.5° W in 
the Gulf of Cadiz, splits into two cores, one at -600 m depth and the other at -1200 m 
depth (not shown), consistent with the observations of Ambar and Howe (1979) and 
Ambar (1982). The two cores follow the northern coast of the Gulf of Cadiz. When the 
cores reach Cabo de Sao Vincente, they turn poleward. A cross-section shows that the 
two deeper cores at 600 and 1200 nn depth (Figure 18a) have core velocities of -4-7 cm/s 
and are relatively salty (Figure 18b) compared to the fresher water above and below the 
cores (The fresher water with depth is due to the general decrease of salinity with depth 
in the region). 

At -38° N a third shallower poleward core is also seen with a core velocity of -4 
cm/s at -400 m depth (Figure 18a). Meddies occur as regular features in the vicinity of 
the southwest coast of Portugal. For example. Figure 19a shows Meddies of 0(100 km) 
at -1200 m depth at -39° N, 12° W and at -38° N, 15° W. 

Meddies are also seen to the west of the Gulf of Cadiz. For example. Figure 19b 
shows Meddies of 0(100 km) at -36° N, 11° W, and at -36° N, 15° W. Cross-sections of 
meridional velocity (Figure 20a) and temperature (Figure 20b) at 36° N illustrate the 
anticyclonic rotation pattern for the Meddies. The corresponding cross-section of salinity 
(Figure 20c) shows the relatively salty signature of the Meddies at -1200 m depth. 

In time the Meddies move in a predominantly westward direction with a 
propagation speed of -2 km/day (e.g., compare the position of the Meddies in Figure 19a 
from September with those in Figure 19b from November). The basic characteristics and 
speeds of the Meddies off the southwest coast of Portugal and to the west of the Gulf of 
Cadiz are consistent with the recent observations of Richardson and Tychensky (1998). 


3. Comparisons of Model Results with Observations 

Since this experiment incorporated the most realistic features of the NCCS, it is 
useful to qualitatively compare the results of the model simulations with observational 
data. Since this is an idealized process-oriented study and not a model hindcast, we 
cannot make direct comparisons with data; however, we can investigate whether the 


21 



phenomenological model behavior is qualitatively similar to observational data in the 
NCCS. Overall, the results of this experiment’s complex flow regime highlight the major 
characteristics of the NCCS with relatively close similarities to field observations. 

a. Comparison of Ocean Currents 

The major characteristics of the modeled currents off Portugal and Morocco are 
consistent with the few available field observations and studies. The model develops a 
mean, equatorward, yearlong surface flow (i.e., the CC) of ~ 10-30 cm/s throughout the 
model domain (e.g., Figure 16), consistent with Wooster and Reid (1963) and 
observations by DMAH/TC (1988). A surface equatorward jet with a mean velocity of 
-20-35 cm/s also develops near the coast. 

Below -200 m depth near the coast, the model establishes a poleward 
undercurrent, consistent with Meincke et al. (1979). Off Iberia, the poleward 
undercurrent in Experiment 4 was shown to transport salty MO between -500 and 1300 
m depth (e.g.. Figure 15c), consistent with Fiuza (1982). 

In the experiments with extensive coastline irregularities i.e.. Experiments 3 and 
4, multiple cores (as many as three) in the undercurrent were generated, particularly off 
the southwestern coast of Iberia, consistent with the observations of Ambar and Howe 
(1979) and Ambar (1982). Since there were no multiple cores simulated in Experiments 
1 and 2, this suggests that the coastline irregularities can play a role in the development of 
multiple cores. The simulated areas of the capes off southwestern Iberia were areas of 
enhanced upwelling (e.g., see Figure 10a for the west and south coasts). Since capes are 
areas where the alongshore component of the wind stress is at a local maximum, these 
areas also have maximum surface coastal current velocities (Batteen, 1997). From a 
continuity perspective, it is argued that enhanced subsurface current velocities should also 
be present in these regions. These subsurface velocities are seen as the multiple core 
signature in the undercurrent. 

Consistent with Frouin et al. (1990) and Haynes and Barton (1990), our winter 
results show a coastal surface poleward current (i.e., the Iberian Current) -25-40 km 
wide, which reaches speeds of -17.5 cm/s north of 41° N (e.g.. Figure 17). 
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Our winter results also suggest the presence of a relatively weak coastal 
equatorward undercurrent in the poleward part of the model domain (e.g.. Figure 17). 
While there are no available observations in the NCCS to confirm this feature, it should 
be noted that an equatorward undercurrrent (and also a coastal poleward surface current) 
with similar speeds has been observed by Hickey (1979, 1998) in the fall and winter off 
the coast of Washington, the poleward end of the California Current system. 

b. Comparison of Eddies 

The model’s generation of numerous eddies on the order of 100 km throughout 
the NCCS is consistent with high resolution satellite imagery and observations by Fiuza 
(1984). Preferred locations of eddy generation occur offshore of prominent changes in 
coastline orientation, i.e., off Cabo da Roca, Cabo de Sao Vincent, the Gulf of Cadiz, and 
Cape Beddouzza and Cape Ghir. A comparison of the results of Fiuza (1984) with the 
results of the model simulations shows that the model generates eddies in similar 
locations. 

The salty characteristics and structure of the anticyclonic eddies seen at -1200 m 
depth in the model simulations (e.g.. Figure 20) are consistent with observations of 
Meddies in the Canary Basin (e.g., Richardson et ah, 1989). The locations, propagation 
speeds, and propagation directions of the Meddies (e.g.. Figure 19) are also consistent 
with the recent synthesis of observations in the Canary Basin of Meddies by Richardson 
and Tychensky (1998). 

Like the previous experiments , the preferred locations for Meddy generation are 
in the vicinity of prominent changes in coastline orientation, i.e., off southwestern 
Portugal for the Meddies. It is suggested that the Meddies are generated by the same 
mechanism as in the previous experiments, i.e., due to the basic instability of the 
equatorward coastal jet and the poleward undercurrent. The major difference between the 
Meddies and the other eddies is that the Meddies are generated at depth rather than near 
the surface. This is likely due to the relatively salty MO core of the undercurrent below 
-500 m depth which makes the water in the core saltier than the water above and below 
it. This could lead to a local baroclinic instability process, resulting in the production of 
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Meddies. This scenario is consistent with the results of Kase et al., (1989), who observed 
eddies in the MO off Iberia and postulated that baroclinic instabilities generated by the 
undercurrent could lead to the generation of Meddies. 

c. Comparison of Upwelling and Filaments 

Upwelling varies by season throughout the domain as a result of the seasonally 
oscillating winds seen in Figure 2. Observations show that upwelling appears first off 
Cape Ghir and Cape Beddouzza in early spring as seen in Van Camp et al. (1991), and 
then progressively to the north as the Azores High moves up the Moroccan coast. This 
progression of upwelling is well simulated in the model. Observed upwelling off the 
west coast of Portugal and along the southern coast of Iberia (Folkard et al., 1997) is also 
well simulated by the model (e.g.. Figure 16). 

Several studies by Fiuza (1983), Barton (1986), and Fiuza and Sousa (1989) have 
confirmed the existence of cold filaments extending off the west coast of Portugal which 
is depicted by the model (e.g., off Cabo da Roca in Figure 16). While investigating the 
waters of southern Iberia in the Gulf of Cadiz, Folkard et al. (1997) observed cold, 
upwelled filaments extending southward from the coast. The model results simulate these 
filaments as well (e.g., south of Cabo de Sao Vincente in Figure 16). 

To the south, satellite observations from Van Camp et al. (1991) and Hagen et al. 
(1996) show a filament of cold, upwelled water extending off Cape Ghir for -200 km. 
They also show that the coastal jet and subsequent frontal zone caused by the adveeted 
cold water can result in the generation of a cyclonic (anticyclonic) meander or eddy on the 
southern (northern) side of the jet. Consistent with these observations, the model depicts 
a filament off Cape Ghir as an equatorward meandering coastal jet which transports cold, 
upwelled water offshore (Figure 16). The model results also show that, consistent with 
the observations, the jet at -32.5° N has a cyclonic eddy (anticyclonic meander) to the 
south (north). 
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IV. SUMMARY 


The objective of this process-oriented study was to use a multilevel PE model to 
investigate the roles of classical as well as unique features in the predominantly wind- 
driven NCCS. Toward this end, results from four numerical experiments of increasing 
complexity were examined. Similar to a classical EEC system, all of the experiments 
produced an offshore surface equatorward meandering current, a surface equatorward 
coastal jet, a surface poleward coastal current, a poleward undercurrent, upwelling, 
meanders, eddies, and filaments. 

To isolate the response to seasonal wind forcing from coastline variations, 
Experiment 1 used seasonal wind forcing along a straight coastline. The results, as 
expected in a classical EEC system, showed that wind forcing was the key generative 
mechanism for the currents and upwelling as well as for the intense and complex 
meander, eddy, and filament structures in the NCCS. 

To examine the additional effects of irregular coastline geometry, two additional 
experiments were run. Experiment 2 used an idealized “straightened”coastline with two 
prominent capes and an embayment representing the Gulf of Cadiz. Due to the presence 
of the Gulf of Cadiz, large anticyclonic meander and eddy formation in the Gulf of Cadiz 
were regular features. Due to the presence of the two prominent idealized capes (i.e., 
Cabo de Sao Vincente and Cape Ghir), cyclonic meanders and eddies tended to form in 
the vicinity of the capes. 

The additional effects of using more realistic capes with extensive coastline 
irregularities in Experiment 3 showed that irregularities in the coastline geometry were 
key elements for enhanced growth of meanders and eddies, and for “anchoring” 
upwelling and filaments to geographical features. For example, upwelling was present in 
Experiment 3 south of Cabo de Sao Vincente and off Cabo de Santa Maria, which was 
not seen in Experiment 2 with the idealized east-west coastline orientation. Experiment 3 
also showed that multiple cores of the undercurrent could be present in the vicinity of 
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capes, e.g., between Cape Ghir and Cape Beddouzza, and between Cabo de Sao Vincente 
and Cabo da Roca. 

In Experiment 4, the additional effects of thermohaline gradients and the MO 
were used. Due to thermohaline gradients, there was an onshore geostrophic flow 
resulting in a tighter temperature gradient all along the coast, and the seasonal 
development of an equatorward undercurrent in the poleward end of the model domain. 

The effects of the MO were seen between -500 to 1300 m depth as relatively salty 
cores of the undercurrent which, west of -7° W, followed the south and west coasts of 
Iberia. At these depths. Meddles occurred as regular features in the vicinity of the west 
coast of the Gulf of Cadiz and off the southwest coast of Portugal. It was proposed that 
baroclinic instabilities of the undercurrent, which had saltier MO cores than the water 
above or below, was the generation mechanism for the Meddles. 

A comparison of the currents, upwelling, meanders, eddies, jets, and filaments 
generated by the model with available observations showed that the model successfully 
simulated the location, size, and velocity of these features. Of the features generated, the 
cores of the poleward undercurrent are of particular interest. Until now, it has generally 
been accepted that multiple cores were the result of MO influence (e.g., Ambar, 1982). 
This study showed, however, that irregular coastline geometry, particularly capes, could 
play an important role in the generation of the subsurface cores. 

In summary, the results from the complex flow regime of Experiment 4, which 
has the most realistic features, highlight the major characteristics and unique features of 
the NCCS with relatively close similarities to field observations. These results show that 
the model successfully captures the qualitative nature of the NCCS. 
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Figure 1. The model domain for the northern Canary Current System (NCCS) is 
bounded by 30° N to 42.5° N, 5° W to 17.5° W. The cross-shore (alongshore) resolution 
is 9 km (11 km). Geographic locations and prominent features are labeled. 
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30m Salinity Forcing at 17.0W: North(42.5N), South (30.0N) 
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Figure 3. Times series plots of monthly temperature and salinity fields used as seasonal 
forcing in the basic simulation. The symbol represents data at 30° N, 17.5° W, while 
the ‘+’symbol represents data at 42.5° N, 17.5° W for 30 m depth. 
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Figure 4. Temperature contours and velocity vectors for Experiment 1 at 30 m depth at 
days (a) 210, (b) 270, (c) 330 and (d) 621. The contour interval is 1°C. To avoid clutter, 
the velocity vectors are plotted every third grid point in the cross-shore direction and 
every fourth grid point in the along-shore direction. Maximum current velocity is 100 
cm/s. 
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Figure 5. Cross-shore sections of meridional velocity (v) for Experiment 1 on day 210 at 
35.5 ° N. Solid lines indicate poleward flow, while dashed lines indicate equatorward 
flow. The contour interval is 2.5 cm/s for poleward flow and 5 cm/s for equatorward 
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LATITUDE : 38N 



LO N G i I U D L 


Figure 7. Time-averaged plots (Experiment 1) for the peak upwelling season (July- 
September) of a cross-shore section of mean meridional velocity(v) in (a) at 38° N, and 
velocity vectors (b) at 1226 m. In (a) solid lines indicate poleward flow, while dashed 
lines indicate equatorward flow. The contour interval for poleward flow is 2.5 cm/s. 
The contour interval for equatorward flow is 5 cm/s in (a). The maximum current velocity 
is 20 cm/s in (b). 
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Figure 10. Temperature contours and velocity vectors for Experiment 3 at 30 m depth at 
days (a) 210 and (b) 285. The contour interval is 1°C. Maximum current velocity is 100 
cm/s. The abbreviations used to denote the different capes are as follows: CR (Cabo da 
Roca), CS (Cabo de Sines), CSV (Cabo de Sao Vincente), CSM (Cabo de Santa Maria), 
CB (Cape Beddouzza), CG (Cape Ghir). 
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Figure 11. Cross-shore sections of meridional velocity (v) for Experiment 3 on day 240 
at 39.5° N (a) and at 35.5° N (b). Solid lines indicate poleward flow, while dashed lines 
indicate equatorward flow. The contour interval is 2.5 cm/s for poleward flow and 5 cm/s 
for equatorward flow. 
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Figure 12. Time-averaged plot of temperature contours and velocity vectors for 
Experiment 3 at 30 m in the third year for August. The contour interval is 1° C. The 
maximum current velocity is 100 cm/s. 
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LONGITUDE 


Figure 13. Time-averaged plots (Experiment 3) for cross-shore sections of mean 
meridional velocity in August at 37.5° N (a) and in January at 39.4° N (b). Solid lines 
indicate poleward flow, while dashed lines indicate equatorward flow. The contour 
interval is 2.5 cm/s for poleward flow and 5 cm/s for equatorward flow. 


53 
















LATITUDE 



Figure 14. Temperature contours (a,b) and velocity vectors (a-c) for Experiment 4 at 30 
depth on day (a) 30, (b) 255, (c) 363. The contour interval is 1° C. Maximum velocity is 
100 cm/s in (a,b) and 20 cm/s in (c). 
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Figure 15. Cross-shore sections of meridional velocity (v) (a-b) and salinity (S) (c) for 
Experiment 4 on day 210 at 37.5° N (a) and on day 363 at 41.5° N (b-c). The contour 
interval for v is 2.5 cno/s for the poleward flow in (a,b), and 5 cm/s for equatorward flow 
in (a) and 2 cm/s in (b). The contour interval for S is 0.01 (0.05) above (below) 1400 m 
depth. 
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Figure 16. Time-averaged plot of temperature contours and velocity vectors for 
Experiment 4 at 30 m in the third year for August. The contour interval is 1° C. The 
maximum current velocity is 100 cm/s. 
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Figure 20. Cross-shore sections for Experiment 4 on day 1059 (end of November) at 36° 
N of meridional velocity (v) (a), temperature (T) (b), and and salinity (S) (c). The contour 
interval for v is 2.5 cm/s for poleward flow and 5 cm/s for equatorward flow. The 
contour interval for T is 0.2 (1.0) above (below) 600 m depth. The contour interval for S 
is 0.01 (0.05) above (below) 1300 m depth. 
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Table 1. Summary of specific experimental design. 


Cases 

Type of coastline 

Mediterranean Outflow 
and 

Thermohaline Forcing 

1 

Straight 

NO 

2 

Ideal Irregular 

NO 

3 

Realistic Irregular 

NO 

4 

Realistic Irregular 

YES 


71 




















Table 2. Values of Constants Used in the Model 


Constant 

Value 

Definition 

~To 

278.2°K 

Constant Reference Temperature 

^0 

34.7 

Constant Reference Salinity 

Po 


Density of Sea Water At Tq and Sq 

a 

2.4 X 10-^ (°K)"^ 

Thermal Expansion Coefficient 

P 

7.5 X lO"* 

Saline Expansion Coefficient 

K 

10 

Number of Levels In Vertical 

Ax 


Cross-Shore Grid Spacing 

Ay 

1.1 X 10^ cm 

Alongshore Grid Spacing 

H 


Total Ocean Depth 

A/ 

800 s 

Time Step 

fo 


Mean Coriolis Parameter 

8 


Acceleration of Gravity 


2 X 10*’ cm'* s’* 

Biharmonic Momentum Diffusion Coefficient 

Ah 

2 X 10*’ cm'* s’* 

Biharmonic Heat Diffusion Coefficient 

Km 

0.5 cm^ s’* 

Vertical Eddy Viscosity 

Kh 

0.5 cm’ s’* 

Vertical Eddy Conductivity 
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Table 3: Salinity profiles 


Layer 

(Depth in meters) 

Salinity (psu) 

1 

10 

36.07 

2 

30 

36.05 

3 

75 

36.05 

4 

150 

36.05 


250 

36.05 

6 

400 

35.75 

7 

600 

35.75 

8 

1226 

36.25 

9 

2283 

35.25 

10 

3656 

34.70 
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